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ABSTRACT 

The first set of supermassive black hole mass estimates, published from 1977 to 1984 by E. A. Dibai, 
are shown to be in excellent agreement with recent reverberation-mapping estimates. Comparison 
of the masses of 17 AGNs covering a mass range from about 10 6 to 10 9 Af Q shows that the Dibai 
mass estimates agree with reverberation-mapping mass estimates to significantly better than ±0.3 
dex and were, on average, only 0.14 dex (~ 40%) systematically lower than masses obtained from 
reverberation mapping. This surprising agreement with the results of over a quarter of a century ago 
has important implication for the structure and kinematics of AGNs and implies that type-1 AGNs 
are very similar. Our results give strong support to the use of the single-epoch-spectrum (Dibai) 
method for investigating the co-evolution of supermassive black holes and their host galaxies. 
Subject headings: galaxies:active - - galaxies: quasars: general - ■ black holes:masses - - galax- 
ies : active: variability 



1. INTRODUCTION 

It is exactly 100 years this year since the public ation 
of th e first evidence of nuclear activity in galaxies (jFathl 
1908). Over the last half century or so, active galactic 
nuclei (AGNs) have been the subject of increasingly in- 
tensive study which has resulted in tens o f tho usands of 
paper s being published. IZel'dovid] (|1964h and ISalpeterl 
(1964) proposed that the huge energy release lasting for 
millions of years from a typical AGN could be explained 
by the accretion of matter onto a supermassive black 
hole. In such accretion, the energy output efficiency may 
reach 43% of the accreting matter's rest mass energy. 
This gave the first estimates of lower limits to the masses 
of AGNs because the luminosity, L, cannot go much 
above the Eddington limit, L E dd ~ 1.3 x 10 38 (A//M o ) 
erg/s, so the mass of the black hole in an AGN has to 
be of the order of several million to several billion so- 
lar masses, depending on t he luminosity of the AGN 
(jZerdovich fe Noviko 
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While the Eddington limit gives a lower limit to the 
mass, Mbh, of the black hole, the actual mass could be 
orders of magnitude greater than this. To understand 
the working of AGNs, Mbh needs to be determined ob- 
servationally, so estimating Mbh has always been con- 
sidered to be a matter of utmost importance in AGN 
studies. If the motions of gas clouds are dominated by 
gravity, masses can be estimated in principle from the 
virial theorem if we know a v elocity and an appropriate 
distance from the center (e.g.. lWoltjerill959l ). The veloc- 
ity along the line of sight can easily be determined from 
the Doppler broadening of emission lines, but determin- 
ing the distance of the emitting material from the central 
object is difficult. 

It was not until the work of lDibail (|1977l 119781 11980L 
1981, 1984a, b) that an attempt was made at a consis- 
tent spectroscopic determination of the masses of the 



central objects for a large sample of AGNs. This en- 
abled Dibai to determine black hole masses and Edding- 
ton r atios, L/LE d d, for dozens of AGNs for the first 
time dDibail Il980. 19843), and to be gin investigations 
in what promised to be (and indeed has turned out to 
be) a very fruitful area: the relationships between these 
quantities an d ot her AGN parameters (see, for example, 
iDibail Il984al and iDibai fc Zasovj 119851) . Unfortunately, 
this work was cut short by Dibai's premature death, and 
in the two decades after his 1977 paper there were only a 
few papers by others us ing the Dibai method to estimate 
black hole masses (e.g., jjolf et al. 19851 iWandel fc Yahill 
[1981 [Padovani fc Rafanellilll988D . 

In the last decade, however, there has been an enor- 
mous growth of interest in determining masses by Dibai's 
method because of the close relationship between the 
masses of black holes and the masses of the bulges of 
their host galaxies (see iKormendv fc Geb hardt 200ll for 
a review). It is only through the Dibai method that 
large numbers of black hole masses in AGNs can cur- 
rently be determined, and the method has already been 
used in making tens of thousands of mass estimates for 
AGNs of all redshifts in the Sloan Digital Sky Survey 
/e.g.. iMcClure fc DunlorJ l200i iSalviander et alJ [20071 : 
iGreene fc Hd 120071: iShen et al.l 12008ft . It is therefore of 
interest to revisit the original Dibai mass estimates and 
see how they compare with more recent independent es- 
timates; 

Diba i fc Pronikl (|1967ft showed that the emitting re- 
gions of the broad and narrow components of AGN 
optical emission lines (what we now call the BLR 
and NLR respectively) had to originate in different 
locations in space, and that the emitting gas, in 
both cases, had to have a cloudy structure, that is, 
to fill only a small fraction, e, of the volume of 
the corresponding region. It has long been recog- 



2 



BOCHKAREV AND GASKELL 



nize d (e.g.. iBahcall. Kozlovskv fc Salpeterl Il972t see 
also IBochkarev fc Pudenkol 119751 ) that BLR variability 
timescales could be used to estimate the effective dis- 
tance of the emitting gas from the black hole, and hence 
t o get masses via the virial theorem. 

iPronik fe Chuvaevl (|1972h pr esented the first long-term 
H/3 li ght curv e for an AGN. iLvutvi fc Cherepashchukl 
(| 19721 ) and ICherepashchuk fc Lvutvil (|1973l ) made 
narrow-band observations of three AGNs over 
several months and published the first actual 
estimates of BLR sizes from the time lag be- 
tween continuum va ria bility a nd line variability. 
IBochkarev fc Antokhinl (11982ft. iBlandford fc McKed 



(I1982D. ICapriotti. Foltz. fc Peterson] (|1982D . and 
lAntokhin fc Bochkarevl (|1983f) independently developed 
methods for recovering information about the BLR from 
the response of the lines to continuum variations, a 
subject which has now become known as "reverberation 
mapping" . Determining BLR sizes became practical and 
widespread with the introduction of cross-correlation 
techn iques (|Gaskell fc Sparkel 119861 ; iGaskell fc Peterson] 
1987) a few years later. 

Once BLR radii were being reliably obtained from 
cross-correlation studies, the main remaining problem in 
estimating masses was in establishing that the gas mo- 
tions were dominated by gravity. The dominant belie f 
from the early days of AGN studies (e.g.,[BurbidgejT958|) 
was that emission- line gas was outflowing from AGNs, 
and the virial theorem obviously cannot be applied to an 
outflowing wind. The idea of determining BLR kinemat- 
ics diagnostics by line-pro fi le var iations was first brought 
forward by S. N. iFabrikal {1980), then Dibai's graduate 
stude nt. The first velocity-resolved reverberat ion map- 
ping (lGaskeTilll988t iKoratkar fc Gaskelll Il989t see also 
IShapovalova et alj l2001al lbl) showed that the BLR was 
not outflowing, but instead showing some net inflow 
in c ombination with Keplerian and/or chaotic motion 
(see IGaskell fc Goosmannll2008D . This discovery imme- 
diately permitted the first rever beration-map ping deter- 
minations of black hole masses HGaskellll988h . 

Despite the promis ing results of the pi oneer - 
ing observations of |Lvutvi fc Cherepashchukl (|1972D 
and ICherepashchuk fc Lvutvil (|1973l ).~ it was em- 
phasized by Bochkarcv ( 198 j), confi r ming earlier 
calculations of IBochkarev fc A ntokhin ( 1982) and 
lAntokhin fc Bochkarevl (|1983f) . that reliable results could 
only be obtained from of long time-sequences of well- 
sampled, high-accuracy, spectral and photometric ob- 
servational data obtained through a large-scale interna- 
tional project. The urgent need fo r starting such coll ab- 
orations was discussed in detail by IBochkarev! (|1987al fbl). 

The following decades saw the progress of the Inter- 
national AGN Watch (I A W) program (see IClavel et al.l 
I199U IPeterson et al.lll99ll et seq.) aimed at determin- 
ing the BLR size and structure from reverberation map- 
ping. That project proved to be one of the largest global 
astronomical monitoring programs to date. More than 
200 astronomers from 35 countries cooperated for 15 
years in accumulating long, densely-sampled, UV, opti- 
cal, and other wavelength time series for many AGN. 
As a result of this, and of additional optical moni- 
toring campaigns, reverberation-mapping estimates of 
the AGN central black hole mass have now been ob- 
tained for over 40 AGNs (see IPeterson et all [2004 and 



IVestergaard fc Petersonll2006f) . 

In this paper we make a comparison of AGN central ob- 
ject mass values from reverberation mapping campaigns 
with the masses obtained by Dibai over a quarter of a 
century ago. In Section 2 we briefly describe the as- 
sumptions made by Dibai. In Section 3 we compare the 
results yielded by the two methods. Section 4 discusses 
the implications of the comparison. 

2. THE DIBAI METHOD OF ESTIMATING MASSES 

In his first (1977) paper Dibai explained the basic prin- 
ciples of his method for estimating the masses of the 
black holes in AGNs and gave masses for 15 AGNs. In 
1980 and 1 981 further masse s for many more AGNs were 
published ( Dibai [T9801 fl981). He then composed a cat- 
alogue of the main characteristics of 77 Seyfert 1 nearby 
galaxies and nearby quasars and of 24 Seyfert 2 galaxies, 
and used it as the basis for new AGN mass and accretion 
rate estimat es. What wer e to be the final results were 
published in Dibai (1984a.b), which did not appear until 
after the author's untimely death occurred in November 
of 1983. 1 

Dibai estimated the mass of the central region of AGNs 
under the assumption that the gas clouds responsible for 
broad emission line formation were moving with more or 
less parabolic velocities in the gravitational-field of the 
black hole so that 



Mbh = l-5Rv 2 /G. 



(1) 



Here Mbh is the mass of the central object (presumably 
a black hole), R is the BLR radius, v is the gas velocity 
(which Dibai determined from the FWHM of the broad 
component of H/3), and G is the gravitational constant. 

To estimate the size of region producing the broad com- 
ponent of H/3 Dibai used the simple relationship: 



e(47ri? 3 /3) = L Hft /E(n,T), 



(2) 



where e is the BLR gas filling factor; that is, the frac- 
tion of the BLR volume filed with the emitting gas. It 
should be n oted that s ince L ex Lhb over many orders of 
magnitude (|Yedll980t IShudenll98lj ). Eq(2) means that 



tfoci 1 / 3 . 



(3) 



Dibai estimated the volume occupied by the emitting 
gas by assuming that the H/3 line is optically thin and 
emitted in the low-density approximation (the so-called 
"coronal approximation"); that is, that the results ob- 
tained for classical ionized hydrogen zones may be ap- 
plied to the line. 

On the right-hand side of Eq. (2), L{Hj3) is the lumi- 
nosity in the H/3 line and E(n, T) = 1.21 x 10~ 7 erg cm" 3 
s" 1 is emissivity for H II zones heated up to T — 10 4 K 
with electron number density, n e = 10 9 cm~ 3 . Dibai 
based his estimate of the density on two things. Firstly, 
the absence of broad components of forbidden spectral 
lines in AGN spectra leads to a lower limit of n e > 10 s 
cm~ 3 . Secondly, the presence of the semi- forbidden CIII] 

1 It should be mentioned that significant parts of the spectral 
and photometric data used by Dibai in above the mentioned papers 
were obt ained with Dibai's personal par t icipat ion in observations: 
e.g., see IDibai. D oroshcnko. & Postno^ l|1981h for results of pho- 
tometry for 27 AGNs and the series of 8 papers by Arakelian, Dibai 
and Esipov with AGN spectroscopy published in t he journal As - 
trophysics in the period 1970-1973 (see references in IDibai l l 1984bD . 
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TABLE 1 

COMPARISON OF MASSES ESTIMATED BY DIBAI 
WITH MASSES ESTIMATED BY REVERBERATION 
MAPPING 



Name 


Dibai 


Reverb 




± 


Ref. 


Dibai-Re\ 


Mrk 335 


7 


50 


7 


15 





12 


2 


0.35 


PG 0026+129 


8 


15 


8 


59 





11 


1 


-0.44 


F 9 


8 


00 


8 


11 





10 


2 


-0.41 


Mrk 590 


7 


10 


7 


68 





07 


2 


-0.28 


3C 120 


8 


00 


7 


74 





21 


2 


0.26 


Ark 120 


8 


00 


8 


18 





06 


2 


-0.18 


Mrk 79 


7 


50 


7 


72 





12 


1 


-0.22 


Mrk 110 


7 


20 


7 


40 





11 


2 


-0.20 


NGC 3516 


7 


40 


7 


63 





15 


1 


-0.23 


NGC 3783 


7 


10 


7 


47 





08 


2 


-0.37 


NGC 4051 


6 


00 


6 


28 





19 


1 


-0.28 


NGC 4151 


7 


20 


7 


12 





16 


2 


0.08 


3C 273 


8 


50 


8 


95 





09 


2 


-0.45 


Mrk 279 


7 


90 


7 


54 





12 


2 


0.36 


NGC 5548 


7 


70 


7 


83 





02 


2 


-0.13 


Mrk 509 


7 


70 


8 


16 





01 


2 


-0.45 


NGC 7469 


7 


30 


7 


09 





05 


2 


0.21 



1 Peterson et al. (2004) 

2 Vestergaard & Peterson (2006) 



9.5 




5.5 



5.5 6.0 6.5 7.0 7.5 8.0 8.5 9.0 9.5 

log M Dibai 

Fig. 1. — Comparison of the mass estimates of IDibail Jl984b), 
M Dibai with reverberation mapping masses, M rever i,. The solid 
line is the OLS-bisector fit and the two dotted lines show regres- 
sions of (X— Y) and (Y— X). 

A1909 line implies an upper limit of n e < 10 10 cm~ 3 . For 
the filling factor parameter, e, Dibai adopted a value ob- 
tained for the fraction of the volume of the emitting gas 
in the Crab Nebula, e = 10~ 3 . With these assumptions 
Dibai then estimated Mbh f° r almost 80 type-1 Seyfert 
galaxies and nearby quasars. 



3. BLACK HOLE MASSES IN AGN: DIBAI (1984B) VS. 
REVERBERATION MAPPING 

Masses derived from more than 15-ye ars of reverber- 
ation map ping of 35 AGN are given in IPeterson et al.l 
(j2004h and lVestergaard fc Peterson! (|2006l ) . 17 AGNs ap- 
pear in both these lists and Dibai's lists. The masses 
range from M BH ~ 1O 6 M up to ~ 10 9 M Q . The black 
hole mass estimates for the AGNs in common are shown 
in Table 1 and plotted in Fig. 1. 

As can be seen from Fig. 1, there is obviously a good 



correlation. Since there are errors in both axes we have 
performed a n ordinary least s quares (OLS) bisector re- 
gression fsee llsobe et aI1ll990D . This gives the relation- 
ship: 

\ogM reverb = l.UlogM Dibai - 0.85. (4) 

As can be seen from Fig. 1, there is obviously a good 
correlation. Monte Carlo simulations show that the error 
in the slope of the OLS-bisector line is ±0.18, so the slope 
does not differ significantly from unity. The systematic 
difference in the logarithms of the masses determined by 
the two methods is 

< \ogM reverb - log M D ibai > = -0.14 ± 0.07. (5) 

The systematic difference between the two methods is 
therefore only 40%. 

The scatter in the ratio of masses determined by the 
two methods is ±0.28 dex, but at least some of this must 
be due to the errors in the reverberation-mapping es- 
timates. The mean of the f ormal reverberation mass 
measurement errors quoted bvlPeterson et all (|2004fl and 
IVestergaard fc Peterson] (|2006l ) is ±0.10 dex. If we make 
the (unlikely) assumption that these measurement errors 
are the only source of error in the reverberation mass es- 
timates, then the average error in the Dibai estimates 
is ±0.26 dex. However, this is an approximate upper 
limit to the errors in the Dibai method because there 
are other (unknown) errors in the reverberation map- 
ping mass determinations. If we assume that the er- 
rors are equally distributed between the Dibai method 
and the reverberation-mapping method, then the mean 
of the errors in the Dib ai method of black h ole mass es- 
timation is ±0.20 dex. iDennev et al.l (|2009[ ) have made 
a detailed examination of the effects of random and sys- 
tematic observational errors on masses estimated by the 
Dibai method. They find that with careful treatment 
of line profiles using high-quality spectra (signal-to-noise 
ratio > 20:1) it is possible to get errors of ^ ±0.10 dex in 
the mass if the observational uncertainties are the only 
source of error. Since Dibai did not have the advantage of 
modern digital spectra, the contribution of observational 
errors to his mass estimates must have been significantly 
greater than ±0.10 dex. Our comparison thus tells us 
that under ideal circumstances, the intrinsic accuracy of 
the Dibai method is potentially quite high. 

4. WHAT DOES THE AGREEMENT TELL US? 

The surprising agreement of the estimated mass val- 
ues implies that all the "classical" type-1 AGNs (Seyfert 
1 galaxies and nearby quasars) are very much alike in 
their properties and structure, and that object-to-object 
differences are smaller than has been hitherto thought. 
The agreement of the Dibai and reverberation mapping 
results provides support for the following: 

1. BLR gas motions are dominated by gravity, as 
has b een shown by velocity-resolved reverber ation map- 
ping (|Gaskelllll98"a iKoratkar k, Gaskel]||1989l . etc.), and 
as is also strongly supported by the inverse correlation 
of line widths with the sizes of the emitting regions of 
different ions, and hence the co nsistency of mass esti- 
mates from these different ions (IKrolik et all HMD: see 
also IPeterson fc Wandell H999. 2000). Thus the emis- 
sion comes predominantly from gravitationally-bound 
gas, and not gas flowing away from the nucleus. This is 
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an assumption in common to both the Dibai and rever- 
beratio n mapping methods. Dibai recognized from the 
outset (|Dibailll977l ) that it had not yet been established 
that the motions were dominated by gravity. He also 
noted that outflow velocities in Wolf-Rayet stars only 
slightly exceed the escape velocity, so his mass estimates 
could also apply if the BLR was a radiatively-driven out- 
flow. 

2. The kinematics of all BLRs are similar. This is 
again an assumption in common to both the Dibai and 
reverberation-mapping methods. 

3. The BLR size scales with optical luminos- 
ity as R oc L a . This is supported by reverbera- 
tion mapping estimates of R (iKoratkar fc Gaskelllll991l : 
IPeterson fc Wandellll999l ; IBentz et al.ll2006f) . Dibai took 
a to be ~ 0.33 (see Eqs. 2 and 3) while IBentz et al.l 
(j2006f ) get 0.52 ± 0.04. While the slope of the line in 
Fig. 1 (1.14 ± 0.18) is already consistent with unity, if 
we adopt a = 1/2, rather than the a = 1/3 Dibai with 
Lh/3 (see Eqs. (2) and (3)), this changes the slope of the 
line to 0.97. 

4. The spectral energy distribution (SED) is very sim- 
ilar in all type-1 AG Ns, as has already been argu ed by 
iGaskell et all (|200l and IGaskell fc Benkeri ((20081) . 

5. The "warm" gas has a similar filling factor, e ~ 
0.001 for all the type-1 AGNs. This means that the space 
between BLR clouds is about 10 times bigger than the 
mean cloud size. Since the BLR is p robably flattened (see 
IGaskell. Klimek. fc Nazarovall2008l ). the average separa- 
tion between the clouds will be a little smaller. 

6. Slow temperature variations of the "warm" (T ~ 
10 4 K ) BLR gas discussed recently by ([Popovic et al.1 
2008) do not make significant deviations of the average 
value) of BLR gas emissivity E(n,T) in H/3 line for the 
17 AGNs considered here from the quantity 1.21 x 10 ~ 7 
erg cm" 3 s _1 adopted by Dibai (1984b). 

7. The possible presence of other BLR emission 
components, such as a BL R component near the jet 
in r ad io-loud AGN (see iBochkarev fc Shapovaloval 
2007t iNazarova, Bochkarev. fc O'Brienl 120071 and 
Arshakian et al.l I2008D in addition to the "standard" 
BLR associated with the accretion disk does not no- 
ticeably affect the average value of the porosity, e, 
of the gas responsible for the BLR. The contribution 
of non-standard BLR components to the total BLR 
emission is thus probably small, at least for the typical 
BLRs of the 17 AGNs considered here. 

The first two of these assumptions are common to both 
the Dibai method and reverberation mapping, and these 
two assumptions tell us things about BLR gas flow in 
AGNs. The next two assumptions are unique to the 
Dibai method and, in combination with the good agree- 
ment of mass estimations we have found, they tell us 
additional significant things about AGNs. The third 
assumption is that the physical conditions in the BLR 
clouds are similar (i.e., similar densities and ionization 
parameters), and the fourth is that the SEDs are simi- 
lar. The fifth assumption was only used to get the scale 
factor or zero point in the original Dibai method. The 
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large number of more recent applications of the Dibai 
method scale the radii to reverberation mapping radii 
instead. Nonetheless, the agreement we find between 
Dibai's mass scale and that of reverberation mapping 
implies that Dibai's assumption of a filling factor of 
e = 0.001 is a good one. The agreement arises, of course, 
because of the agreemen t in estimat ed radii. Dibai in- 
deed noted at the outset (Dibai 1977) that the estimated 
radii agreed with the reverbera tion mapping results of 
iLvutvi fc Cherepashchukl (|1972fl . 

Taken together, the apparent correctness of the main 
assumptions of Dibai's method reinforces the idea that, 
despite the wide range of masses of their black holes and 
the wide range of accretion rates, a large fraction of type- 
1 AGNs are surprisingly similar. This is an important 
result which needs to be understood and which requires 
further study. 

The final, obvious, and important conclusion that 
can be drawn from the good agreement we find be- 
tween masses estimat ed by t he sin gle-epoch spectrum 
method pioneered by iDibail (|1977l ) and by the rever- 
beration mapp i ng m ethod (and also the agreement that 
IPeterson etall 120041 and IVestergaard fc Peterson] 120061 
find between their own single-epoch estimate and rever- 
beration mapping) is that the Dibai method does give 
reliable black hole mass estimates. This is very impor- 
tant for studying cosmic evolution of black holes and 
host galaxies because it shows that the long-term, labour- 
intensive monitoring of reverberation mapping is not nec- 
essary to obtain the mass of every single AGNs in the 
sample. The reliability of the Dibai method makes in- 
vestigations of cosmic AGN evolution, black hole growth, 
and the links with the properties and evolution of host 
galaxies considerably easier. The simple comparison we 
have presented here implies that we can indeed trust this 
method which has now been used already for obtaining 
mass estimates of many tens of thousands of AGNs. 

The success of the Dibai method certainly does not 
mean that there is no need for further reverberation map- 
ping studies. Reverberation mapping is needed for test- 
ing the reliability of black hole mass determinations for 
AGNs of extreme types which are not represented in Ta- 
ble 1. For example: high Eddington ratio AGNs (so 
called "narrow- line Seyfert Is" ) and very low Eddington 
ratio AGNs (such as FR1 radio galaxies), which might 
have non-classical BLRs, are not represented. The good 
agreement of the Dibai method with reverberation map- 
ping for the 17 AGNs considered here does not reduce 
the importance of the reverberation mapping in general. 
The two methods are complimentary and not mutually 
replaceable. 
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